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Possible solution to the 7 Li problem by the long lived stau 
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Modification of standard big-bang nucieosynthesis is considered in the minimal supersymmetric 
standard model to resolve the excessive theoretical prediction of the abundance of primordial lithium 
7. We focus on the stau as a next-lightest superparticle, which is long lived due to its small mass 
difference with the lightest superparticle. It provides a number of additional decay processes of 
7 Li and 7 Be. A particularly important process is the internal conversion in the stau-nucleus bound 
state, which destroys the 7 Li and 7 Be effectively. We show that the modification can lead to a 
prediction consistent with the observed abundance of 7 Li. 
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I. INTRODUCTION 



The theory of big-bang nucleosynthesis (BBN) has 
been successful in predicting the abundance of light ele- 
ments in the Universe from a single parameter, baryon- 
to-photon ratio rj. The recent results of the Wilkinson mi- 
crowave anisotropy probe (WMAP) experiment how- 
ever, put this theory into challenge. The extraordinarily 
precise results from WMAP are put together with the 
standard BBN (SBBN) to predict the abundance of 7 Li 
to be (4.15lg;||) x l(r 10 [l] if we adopt r\ = 6.1 x KT 10 
(68 % C.L.) This prediction is inconsistent with 

the observation of metal-poor halo stars which implies 
(1.23±g:l|) x 10" 10 3 reported by Ryan et al. [|. The 
inconsistency persists even if we adopt the recent ob- 
servations, which give the less restrictive constraint of 
(2.19±£^) x KT 10 1 and (2.34±£^) x 1CT 10 @. This 
discrepancy can be hardly attributed to the correction of 
the cross section of nuclear reaction 0, Q , and astrophys- 
ical solutions are pursued Q. 

Another interesting approach to this problem would 
be to consider the effects induced by new physics be- 
yond the standard model (SM). Exotic particles which 
interact with nuclei will open new channels to produce 
and destroy the nuclei, giving a potential solution to the 
7 Li problem. In this paper we investigate a possibility 
that the interaction is initiated by a formation of the 
bound state of exotic negatively charged massive parti- 
cles (CHAMPs) and a nucleus. (For the other solutions, 
see [H EH E3-) 

So far the bound-state effects by CHAMPs have been 
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bly charged particles, see also Refs [19j. In particular, a 
significant enhancement of a 6 Li-production rate through 
4 Hc + D — > 6 Li + 7 by the bound state with 4 He was re- 
ported [IH for the first time and recently confirmed [2(| • 
This hinders the compatibility between particle physics 
models and BBN [U. 

In addition, some nonstandard effects on the abun- 
dance of 7 Li and 7 Be were also considered in Ref. [l6| 
and more recently in Ref. [1^ . Introducing the CHAMPs 
with the mass of electrowcak scale, the authors in 
Ref. [22| newly considered several destruction channels 
of 7 Be nuclei through the trapping of the CHAMPs to 
show that the abundance of the CHAMPs needs to be 
larger than 0.02 per baryon and that their lifetime has 
an allowed window between 1000 and 2000 sec. 

We put the CHAMP BBN scenario in the minimal su- 
persymmetric standard model (MSSM) with the conser- 
vation of R parity. MSSM doubles the particle content of 
the SM by introducing supcrparticles, which can accom- 
modate the CHAMPs. The CHAMPs need a lifetime 
long enough to sustain the sufficient abundance at the 
time of nucleosynthesis. Although the R parity conser- 
vation stabilizes the lightest superparticles (LSPs), the 
observational constraints exclude charged supcrparticles 
as a candidate for LSPs, which is usually considered to 
be neutralinos x° or gravitinos. A possible candidate of 
CHAMPs is the next-lightest superparticle (NLSP) with 
electric charge, which can have a long lifetime by assum- 
ing a small mass difference from the LSP j23| . 

We assume in the present paper that the LSP is a neu- 
tralino and the NLSP is a stau f, the superpartner of tau 
lepton r. The staus can decay into neutralino LSP with 
the hadronic current, through which they also interact 
with the nuclei. The gravitino LSP, on the other hand, 
does not couple with hadronic current. We consider the 
bound state of 7 Be and f ~ in the early Universe and the 
subsequent decay chain of nucleus 7 Be — ► 7 Li — > ' He due 
to the interactions of the two. The 7 He nuclei rapidly 
decay into 6 He nuclei, which are effectively stable in the 
considered time scale. With the freedom of the mass of 
stau rrif and its lifetime t> , we search for the possible 
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solution to the 7 Li problem that arc phcnomcnologically 
acceptable. 

This paper is organized as follows. In Sec. |TT] we 
overview some new decay channels by stau and estimate 
their lifetimes. Here we will sec that stau-nuclcus bound 
states play an important role in 7 Be/ 7 Li reducing pro- 
cesses. In Sec. IIIII we numerically calculate primordial 
abundances of light elements while taking into account 
the new channels. Then we will see the possible solution 
of the 7 Li problem. Finally, summarization is made in 
Sec. HVl 



II. THE DESTRUCTION OF 7 Be/ 7 Li IN MSSM 

A. Elementary interactions of the staus 

We consider a modification of the SBBN scenario un- 
der the MSSM. MSSM introduces a set of superparticles 
as the partners of the particle appearing in the standard 
model (SM). The superparticles interact with the stan- 
dard particles and thus introduce additional decay chan- 
nels of 7 Be and 7 Li to the standard BBN theory. The 
additional channels give a possible solution to the prob- 
lem where the theoretical prediction of the abundance 
of 7 Be and 7 Li, or collectively 7 Be/ 7 Li, exceeds the ob- 
servational results by a factor of ~ 2-3. We consider 
in the present paper that the destruction of primordial 
7 Be/ 7 Li nuclei is due to their interaction with the neg- 
atively charged staus f ~ , the superpartner of the tau 
lepton t~ , which we identify as the next-lightest super- 
particle (NLSP). The mass eigenstatc of the stau is given 
by the linear combination of the left-handed stau tl and 
the right-handed staus tr. as 



(1) 



where 9 T is the left-right mixing angle and j T is the CP- 
violating phase. 

Staus have attractive features when considering the de- 
struction of 7 Be/ 7 Li. First, staus have a negative charge 
and can form a bound state with nuclei so that they inter- 
act efficiently. Second, staus couples with the hadronic 
current J^, through which they interact with nuclei as 
we see below. Third, staus can be abundant at the time 
of BBN. They can acquire the sufficicntlty long lifetime 
when the staus and LSPs, which we assume as the nu- 
tralinos, have a mass difference tiny enough. 

The interaction of staus is described by the Lagrangian 



AG 

Ant = f *x°(. 9Lj P l + g R P R )r + -zz^PltJ, 
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(I^P l ^){D t ^P l t)+U.c. 



(2) 



where G = 1.166 x 10 5 GeV 2 is the Fermi constant, 
Pl and Pr are the chiral projection operators, I S {e, /i}, 






FIG. 1: Feynmann diagrams of the decay of staus. (a) f 
T X°, (b) r -»• nv T x°, and (c) f -> lviv T x° ■ 



and (?l and g R are the coupling constants. These cou- 
pling constants arc written in terms of the SU(2) L gauge- 
coupling constant g; when the neutralino is bino like, for 
instance, they are given by 



fjL 



3r 
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Sill f^w cos t) T , 



cos 9 W 



sin%smfl T e , 



(3) 
(4) 



where #w is the Weinberg angle. The interaction La- 
grangian © give rise to the following decay channels (see 
Fig. [J): 



f -> llSlV T X°- 



(5) 

(6) 
(7) 



Process (O has a typical lifetime O(10 20 )sec, pro- 
cess © has (10~ 6 -10 2 )sec, and © has (10 2 - 10 12 ) sec. 
Since the BBN takes place (1-100) sec after the big 
bang, the staus will decay entirely before BBN unless 
the channel (O is closed. Our scenario therefore re- 
quires Sm < m T = 1.7 GeV. Note that the chan- 
nel ([6]) also closes when 6m is less than the pion mass 
TO^ ~ 140 MeV. Although the required LSP-NLSP mass 
difference is small compared to the typical mass of LSP 
which is O(100 GeV), it is preferable in attributing the 
dark matter (DM) to the neutralino LSPs since it al- 
lows the LSP-NLSP coannihilation. With this tiny 5m, 
the neutralino naturally becomes a cold dark matter in- 
stead of warm or hot dark matter [24[ even though it 
is produced non thermally. Hence our model is free of 
the constraints from the large-scale structure formation 
of the Universe. 



B. Interactions of staus with 7 Be and 7 Li 

In this section, we consider the stau-nucleus interac- 
tion processes that are relevant to the primordial BBN. 
Three processes are discussed: (1) the hadronic-current 
interaction, (2) stau-catalyzed fusion, and (3) internal 
conversion of stau-nucleus bound state. We consider the 
lifetimes of each process because it is crucial to under- 
stand the impact upon the modification of BBN. 
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5m (GeV) 

FIG. 2: (color online). The lifetime of free stau as the func- 
tions of 5m. Here we take m^o = 300 GeV, 9 T = n/3, and 
7 T = 0. The hadronic decay is dominant for 5m > m n while 
the leptonic decay is exclusively allowed for 5m < m n . 



1. Destruction of nuclei by a hadronic- current interaction 
with free staus 

Staus can interact with the nuclei through the hadronic 
current and thereby alter the BBN processes. The abun- 
dances of 7 Li/ 7 Be are changed by the new decay chan- 
nels: 

r^X° + ^r + 7r ± , (8) 
7T+ + 7 Li 7 Bc, (9) 
7r-+ 7 Be^ 7 Li (10) 
7T- + 7 Li ^ 7 He. (11) 

The process 7r + + 7 Hc — -> 7 Li does not occur since 7 He 
is very unstable, while the pions can be either real or 
virtual; here the virtual pion should actually be regarded 
as a hadronic current propagating between the stau and 
the nucleus. The pions produced in the process © also 
change the proton-neutron ratio and thereby change the 
primordial abundance of the light elements. 

We present the lifetime of the free stau in Fig. [5] as 
functions of 5m [23|. Here we take m^o = 300 GeV, 6 T 

'.'). and j T = 0. 



2. Stau-catalyzed fusion 

Another process to destroy the 7 Li/ 7 Be is nuclear fu- 
sion catalyzed by staus. A nucleus has a Coulomb barrier 
which normally prevents the nuclear fusion, while the 
barrier is weakened when a stau is captured to a state 
bound to the nucleus. The nuclear fusion is thus pro- 
moted by forming a stau-nucleus bound state. The stau 
serves as a catalyst and is left out as the fusion proceeds 
through. 

This stau-catalyzed fusion process provides the follow- 




bound state 



FIG. 3: The Feynmann diagrams of internal conversion of 7 Be 
( 7 Li). 



ing decay channels: 

7 Bc + f -> ( 7 Bcf) + 7, (12) 

7 Li + f -> ( 7 Lif)+7, (13) 

( 7 Bef) + p^ ( 8 Bf)+ 7 , (14) 

( 7 Bcf)+n^ ( 7 Lif)+p, (15) 



( 7 Lif) +p^ f + 2 4 He or -^f + 2D + 4 He. (16) 

The lifetime of the stau-catalyzed fusion is estimated to 
be longer than 1 sec [H[ . We follow Ref. [2l[ to calculate 
the stau-catalyzed fusion rate. 



3. Internal conversion of nuclei in the stau-nucleus bound 

state 

The interaction between a stau and a nucleus proceeds 
more efficiently when they form a bound state (see Fig. 
12]) H2 due to two reasons: (1) the overlap of the wave 
functions of the two becomes large since the stau and par- 
ticle are packed in the small space, (2) the small distance 
between the two allows virtual exchange of the hadronic 
current even if dm < . The stau-nucleus bound state 



decays through the following processes: 

f + 7 Be -> (f 7 Bc) -> x° + v T + 7 Li, (17) 

f + 7 Li -> (f 7 Li) -> x° + v T + 7 He, (18) 

7 Hc^ 6 Hc + n, (19) 



He + background particles — > 3 He, He, etc., (20) 

where the parentheses denote the bound states. We note 
that we introduce not only reaction (|17p . but also reac- 
tion (p~8|) . The 6 Hc nucleus can also decay into 6 Li via 
f3 decay with the lifetime 817 msec. We do not take this 
process into account since this process is much slower 
than the scattering process (|20|) . 

The lifetime of the internal conversion tic is obtained 
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from the Lagrangian ([2]) as 



tic 



(21) 



where \ip\ 2 is the overlap of the wave functions of the 
staus and the nucleus, 



(H = 2Ef \ EBc J dLIPS|(x°^ 7 Li|£ int |f 7 Be)| 

X (27r)V 4) (p f + p Bc - P^o - Pv T - Phi), 



(22) 



and 



dLIPS = J 



(2^) 3 2^ 



(23) 



Here i E {\°,^T, 7 Li} for the process (fT7|) and i £ 
{x°> ^t, 7 He} for (JTHJl . The following approximations are 
applied to evaluate the lifetime further. We estimate the 
overlap of the wave functions in Eq. (|21[) by assuming 
that the bound state is in the S state of a hydrogen like 
atom, and obtain 



(24) 



'mid 



where a nuc i = (1.2 x A 1 / 3 ) fm is the radius of the nucleus. 
The matrix element of the nuclear conversion appearing 
in Eq. (|2"2"|) is evaluated by the ft value of the corre- 
sponding (3 decay obtained from the experiments. The 
experimental ft value is available for 7 Li «-> 7 Be but not 
for 7 Li <-> 7 He, however. We assume that the two pro- 
cesses have the same ft value. As long as we consider 
the quantum numbers of the ground state of Li7 and Hc7 
we can expect a Gamow- Teller transition can take place 
since they are similar with those of He6 and Li6 and we 
know that they make a Gamow- Teller transition. The 
Gamow- Teller transition is superallowed and has a simi- 
lar ft value to the Fermi transition such as 7 Li «-> 7 Be. 

Our new effects have been treated as if 7 Li or 7 Be in its 
bound state would have an effectively new lifetime which 
is caused by the virtual exchange of the hadronic current 
with a stau. Thus this new process is not the two-bodies 
scattering. So, there is no corresponding astrophysical S 
factor in these processes. 

The evaluated lifetimes of reactions (fT7|) and (fT8| un- 
der these approximations are presented in Fig. [4] as func- 
tions of 8m. There we take m^o = 300 GcV, 6 T = 7r/3, 
and 7 T = for both reactions. We find that the life- 
time of the internal conversion process is in the order of 
10 -3 sec. The lifetime of stau- 7 Li bound state diverges 
around Sm = TO7 Li — m7 Bc = 11.7 McV, below which the 
internal conversion is kinematically forbidden. 

As we will see later, the internal conversion processes 
are dominant over the other processes. 




FIG. 4: (color online). The lifetimes of internal conversion 
processes as the function of Sm. Top panel: (f 7 Be) — > x° + 
v T + 'Li, bottom panel: (t 7 Li) — > x° + v T + 7 He. We take 



m^o = 300GeV, 



7r/3, and 7t = in both figures. 



III. NUMERICAL CALCULATION AND 
INTERPRETATION OF THE RESULT 

In this section, we study the effectiveness of new de- 
cay channels on the 7 Li problem by numerical calcu- 
lation. To do this, we choose the abundance of stau 
Yf = rif /s|t=t &ee2ie out at freeze out time and mass dif- 
ference Sm as free parameters, since these values are sen- 
sitive to the abundance of 7 Be and 7 Li. Here s is en- 
tropy density. The number of 7 Li interacted with the 
stau is determined by the number density of the stau. 
The hadronic decay rate of stau is mainly determined by 
the mass difference. The decay rate is also determined by 
the stau mixing angle 9 T , CP violating phase 7 T , and neu- 
tralino mass. As we showed in 23|, however, the effects 
by these parameters are much less than mass difference. 

We estimate the number density of bound states by 
using Saha equation, 



^BS 



/m N ry 3/2 

V 2tt 



= £b/T 



(n N - n BS )(n f - - n BS ). 



(25) 
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6m (GeV) 



FIG. 5: (color online). The constraints from the light-element 
abundance shown in the 5m— Y? plane. The white region is 
the parameter space which is consistent with all the obser- 
vational abundance including 7 Li/H=(1.23t°:25) x iO" 10 St- 
The regions enclosed by dotted (green), dashed (light blue), 
and dash-dotted (purple) lines are excluded by the observa- 
tions on 4 He, D and 6 Li, respectively. The thick-dotted line 
represents a yield value of stau whose daughter particle, neu- 
tralino, accounts for all the dark matter component. Here we 
took rj = 6.1 x 10" 10 , m s o = 300 GeV, T = vr/3 and 7 r = 0. 



Here, ubs, Tif, and n-^ denote the number densities of 
the bound state, the stau, and the nucleus, respectively; 
ton, Eb, and T denote the nucleus mass, binding energy 
of the bound state, and the temperature of the Universe, 
respectively. The Saha equation is valid only when the 
expansion rate of the Universe is much smaller than the 
formation rate of the bound state. This condition is not 
quite satisfied in our case. We will explain that our re- 
sults are qualitatively acceptable at the end of this sec- 
tion. For more detailed discussion, see 0, [2f| . 

In Fig. O wc show the parameter region of dm and Yf 
allowed by the observational light-element abundances, 
where wc take T] = 6.1 x lCT 10 , m^o = 300 GeV, T = tt/3 
and 7 r = 0. The white region is the parameter space 
which is consistent with all the observational abundance 
including 7 Li to hydrogen ratio ( 7 Li/H). The regions en- 
closed by dotted (green), dashed (light blue), and dash- 
dotted (purple) lines are excluded by the observations 
on the mass fraction of 4 Hc (Y p ), deuterium to hydrogen 
ratio (D/H) and 6 Li to 7 Li ratio ( 6 Li/ 7 Li), respectively. 
We adopt the following observational constraints on pri- 



mordial light clement abundances 1 2 

Yp° bs = 0.2516± 0.0040 H2, 
(D/H) obs = (2.82 ± 0.26) x 10" 5 

. obs 



Log 10 ( 7 Li/H) c 
( 6 Li/ 7 Li' 



obs 



= -9.91 ±0.10 0, 
< 0.046 ± 0.022 + 0.84 \34 



(26) 
(27) 
(28) 
(29) 



The thick-dotted line represents a yield value of stau 
whose daughter particle, neutralino, accounts for all the 
dark matter component. This line is given by the yield 
value of dark matter 



Y DM = 3.80 x 10" 



DM' 



v) 



-1 



(30) 



with n DM h 2 = 0.104±0.010 (68% C.L.) , here to DM is 
to^o. Yf must be smaller than this value, in order to pre- 
vent the overclosure of the Universe. The cosmologically 
interesting region is below the line. 

We can find another white region even if we adopt 
the more restrictive value of 7 Li/H shown in Ref. [3]. 
The upper central region is excluded by the observational 
constrains on D/H and 4 He due to charged pions emit- 
ted from decaying staus [3^|. In the current model, no 
hadrodissociation processes of light elements occur. 

The qualitative feature of the allowed region is ex- 
palaincd from the following physical consideration. First, 
staus need to have lifetimes Tf longer than the time re- 
quired to form the bound state of a stau and a Be. 
The required time £f rm, 7 Be is estimated from the bind- 
ing energy i?bin, 7 Bc = 1490 keV as if or m, 7 Bc ~ 10 9 sec • 
(£bin/Bo/kev)-2 ^ O(10 2 sec). Figure [2] shows 5m < 
(100 - 200)MeV for Tf > t {oIm r Be and hence the al- 
lowed region appears only in this region. Second, the 
yield value of staus Yf needs to be large compared with 
that of 7 Li which we denote by Yr u . We estimate Yr u 
from the hydrogen to entropy ratio nn/ s ~ 0(1O _1 °) and 
the constraint (1281) as 



Y, 



7 Li 



obs 



• — ~ 10" 



9.63 x 10 -10 



10" 



(31) 



We thus need Yf > 10 , and again the allowed region 
appears in this region. Third, the excessive destruction 
of 7 Li by the process (fTS)) needs to be avoided due to the 
constraint ([29]) . This condition puts a limit to the forma- 
tion rate of the bound state of a stau and a 7 Li. We then 
need either Tf < tf or m, 7 Li ~ 10 9 sec • (E h ^ 7 T i /k ev)~ 2 ~ 
O(10 3 sec) (here we use E hin r Li = 952keV Hi), or Y f 



Here we have used conservative errors of Y° hB 
discussion in Ref. 12611 . See also the other recent observational 



according to a 



2 About the errors of 6 Li/ 7 Li, see also the discussion in 
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to be small enough. The former condition leads to 
5m > lOOMeV, although this region is subject to the 
strong restriction considered in the previous paragraph. 
The latter imposes an upper limit on Yf in the region 
8m < lOOMeV, and Fig. [5] suggests that this limit is in 
the order of 10~ 20 . The tininess of Yf shows that the 
internal conversion processes (|17|) and (|18p are dominant 
over other processes such as the pion exchange and the 
stau-catalyzed fusion. We can confirm this dominance by 
an explicit calculation of the rates of these processes. 

We used the Saha equation in our calculation although 
the formation rate is comparable to the expansion rate. 
Calculation using the Boltzmann equation will give a 
lower number density of the bound state ties and conse- 
quently shift the allowed region upward in dm < lOOMeV 
in Fig. El 

IV. SUMMARY 

We have investigated a possible solution of the 7 Li 
problem in a framework of MSSM, in which the LSP 
and the NLSP are ncutralino and stau, respectively, and 
have a tiny mass difference of Sm < (100— 200)MeV. The 
staus then survive throughout the BBN era as shown in 
Fig. [21 and provide additional decay processes as men- 
tioned in Sec. Ill Bl to reduce the primordial 7 Li abun- 
dance. 

Taking the three new processes into account, we nu- 
merically calculated the primordial abundance of light 
elements varying the LSP-NLSP mass difference and the 
abundance of stau. Taking 9 T = 7r/3, 7 r = 0, we ob- 
tained the parameter region consistent to the observed 



7 Li abundance. 

Though we have shown that the internal conversion 
is very important for the calculation of the primordial 
abundance, we need to improve our calculation for bet- 
ter accuracy. First, we need to include reaction processes 
such as ( 6 Lif ) — > 6 Hc + x + u T . This process can change 
the prediction of 6 Li and hence change the allowed re- 
gion. Second, we should calculate the number density of 
the bound states not by the Saha equation (j2"S")) but by 
the Boltzmann equation. At the formation temperature, 
the capture rate is less than the expansion rate of the Uni- 
verse. Therefore we will obtain a lower number density 
of the bound states and consequently the upward shift of 
the allowed region. Third, we should explore other val- 
ues of the parameters 9 T , j T , and m-g, which affect the 
lifetime of stau and also those of the bound states. We 
leave these for our future works. 
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